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The structure of the nucleon is studied variationally on the lattice by maximizing the overlap be- 
tween the nucleon and a trial function generated by an interpolating field containing variational 
parameters. We examine the effect of the spatial extent of smeared quark sources, the degree of 
gauge field smearing, the positions of smeared quark sources, inclusion of lower Dirac compo- 
Qh nents and of diquark substructure. Exploratory calculations with quenched Wilson fermions at a 

fH pion mass of 900 MeV achieved overlaps as high as 80%, and there was no evidence of diquark 

substructure. 
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1. Introduction 

Experimental studies of nucleon structure have revealed many features of form factors, parton 
distributions, and generalized form factors to impressive precision, and lattice calculations are 
now beginning to calculate these observables from first principles in agreement with experiment. 
However, a complicated numerical calculation that merely reproduces an experimental number 
does not necessarily, in itself, provide any physical insight as to why and how QCD works to 
produce a given electric or magnetic form factor or a particular decomposition of the nucleon spin. 
Hence, it is desirable to use lattice QCD as a tool to obtain insight into the quark and gluon structure 
of the nucleon. 

Calculation of a nucleon on the lattice occurs in two steps. First, we create a state, which we 
will refer to as a trial state, by acting with an interpolating field with the proper quantum numbers on 
the QCD vacuum. Then, that state is evolved in Euclidean time to project out the physical ground 
state. Since, as described below, it is straightforward to extract the overlap between normalized 
nucleon and trial states, we can use the overlap as a quantitative criterion to systematically improve 
the trial state and thereby obtain insight into key features of the nucleon wave function. 

Since the only practical way to calculate nucleon correlation functions on a lattice is to ex- 
press them in terms of quark correlation function, the choices in sources are the positions of each 
quark, the spatial extent of the quark source achieved by gauge invariant smearing, the degree of 
smoothing of gluon fluctuations achieved by smearing the gauge links included in the source, the 
inclusion of only upper or of both upper and lower Dirac components, and the choice of 7 matrix 
structure used to produce the nucleon quantum numbers. 

Given recent resurgence of interest[l -3] in diquarks [4], it is of interest to look for any ev- 
idence of diquark correlations variationally. We could, in principle, study two possible diquark 
configurations in a nucleon[5, 6], the scalar channel (uCysd), and the vector channel {uCy^d). 
However, at the one gluon exchange level, quarks in the scalar configuration have lower energy 
than quarks in the vector configuration [2], which are thus called "good" and "bad" diquarks re- 
spectively, and we will focus our attention only on the good diquarks and consider sources of the 
form {uCjs d) u. Analogous to the simple picture of a meson comprised of an antiquark and quark 
connected by a flux tube, which leads naturally to Regge trajectories, in the limit in which two 
quarks bound into a point-like diquark in a nucleon, one could also develop a "dog bone" model 
of baryons with a diquark and quark connected by a flux tube. This picture is phenomenologically 
successful(Wilczek). We can explore this physics in our trial function, by allowing for a diquark to 
have a different degree of spatial localization and to be separated from the remaining quark. 

2. Variational Method 

In this study, we used a simple variational approach to study the ground-state wavefunction 
of the nucleon. For a trial source operator /(x), the overlap with the ground state is calculated 
from a fit to the nucleon two-point correlation function [7]. Starting with the correlation function 
in position space: 

c(x,o = (f^|/(x,o AO, 0)1^2) 
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we insert a complete set of states and project onto zero momentum in the usual way to obtain: 

c(0^Ic(x,0 = I^'-^"'|(«|/(o,o)|q)|2^£a„.-^"'. 

X n n 

The energies En and coefficients An are extracted by fitting C{t) to a sum of exponentials, and 
the normalized overlap of our trial source with the nucleon ground state is given by: 

l(o|/(o)|a)|2 Ao ^2.1) 



Lj(«|7(o)|a)p c(o) 

Note that this variational study with a single source differs from the variational approach used 
extensively in spectroscopy [8], which considers a superposition of distinct sources with arbitrary 
coefficients and determines the optimal coefficients by minimizing the energy, rather than maxi- 
mizing the overlap. 

2.1 Correction at / = 

For ^ > 0, the correlator C{t) = {Q\J{t)J{&)\Q) is equal to the transfer matrix element 

C{t) = {Q.\Jt-^' J\Q) . (2.2) 

For / = 0, the Euclidean correlator is no longer equal to the simple matrix element (2.2), but 
rather to the vacuum expectation value of the normal-ordered product of /and / [9], 

Cg(0) = (f2|7V[//]|f2), (2.3) 

where normal ordering is defined by its action on quark operators: in a basis where 
7o = diag(l, 1,-1,-1), 

Aj\- (-\- ^-\^ I ^^r(^)^bp{y) ifa,/3 = l,2 

[-qbpiyjqayw if a, 15 = 3, 4. 

In order to compute the matrix element (2.2) without normal-ordering, we have to sum the Eu- 
clidean expectation value over all contractions. In practice, this is achieved by replacing the naive 
quark propagator at equal time by the one corresponding to a non-normal-ordered expectation 
value[10, 11], 

+ Caia,bjp{^^ny,t') , (2.5) 

where 

- N[qaia{^,t)qj,j^{y,t')] (2.6) 
depends on the fermion action. For the Wilson action [9], 

Caia^bjpi^.t'J.t') = -5ij{l-r4)apB-^\^J)S{t,t') (2.7) 
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Log plot of the two point function data with curve fit. 



comparing diagonal points 



2 -component 
4 -component 



Figure 1: Plot of a typical two-point correlation func- Figure 2: Overlap as a function of quark RMS radius 
tion. (in lattice units). Results are shown for both two- and 

four-component spinors. 



with 



Babi^J) Sab 8 {X J) 

-K^ [c/r^x)5(f+fJ) + C/f-(3;)5(f-fJ) 



=1 " 



(2.8) 



In our calculations with the Wilson action, we include this term to calculate the correct value 
of C(0). However, there is no comparable transfer matrix for domain wall fermions, so we cannot 
presently extend this work to chiral fermions. 



3. This Calculation 

We performed our exploratory calculation on 16^ x 32 quenched Wilson lattices with j8 = 6.0 
and K = 0.1530 (rriji ^ 900M^y) for ensemble sizes of 100 - 200 configurations. The trial sources 
were of the form 

Ja={Up[Cr5]pyDy)Ua. (3.1) 

We varied the number of gauge invariant smearing steps for the quark sources, controlling 
their RMS radius [7] , the number of dirac spinor components (four, or two in the non-relativistic 
limit), the gauge field smearing in the source links, the relative size of quark and diquark smearing, 
and the relative position of the quark and diquark. 

3.1 Results 
Quark Smearing 

We smear all three quark fields to the same RMS radius and locate them at the same spatial 
position. Results for the variation of the quark RMS radius are shown in Fig. 2. The overlap 
behaves as in previous studies, starting at the order of 10~^ for a small quark radius, increasing to a 
maximum at some finite radius, and falling off at larger radii. This reflects the finite spatial extent 
of the quark wavefunction within a nucleon. The peak occurs around 4.5 lattice units (^ 0.4 fm), 
consistent with [7]. 
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Figure 3: Overlap as a function of quark RMS ra- 
dius for four-component spinors. Results are shown 
for calculations with and without APE smearing. 



four components, peaks at .36 




Figure 4: Overlap as a function of quark RMS radius 
for two-component spinors. Results are shown for cal- 
culations with and without APE smearing. 



Upper components only, peaks at .51 




Figure 5: Contour plot of the overlap as a func- 
tion of both quark and diquark RMS radius, for four- 
component spinors. 



Figure 6: Contour plot of the overlap as a func- 
tion of both quark and diquark RMS radius, for two- 
component spinors. 



Two vs Four Spinor Components 

When constructing the nucleon two-point function, we can use either the full quark spinor q, 
or the projected spinor which corresponds to the upper two components in the Dirac basis. 

We observe (Fig. 2) that using the two-component spinors increases the overlap significantly. This 
is consistent with the expectation that the lower components be in a p-wave type state, which has 
very poor overlap with the approximately gaussian wavefunction used as a trial source (also see 
discussion in [12]). 

Gauge Field Smearing 

Because of the large gluon fluctuations associated with the link variables, the overlap is in- 
creased by smearing the gauge links used in constructing the quark sources. For this, we used 25 
iterations of APE smearing with e = 0.35 in the notation of [13]. Figs. 3 and 4 show the overlap 
calculated with APE smeared links compared with the overlap for no link smearing. We note that 
smearing the gauge fields significantly improves the calculated overlaps. 

Relative Size of Quark, Diquark 

Using different fermion smearing parameters for the diquark and lone quark parts of the nu- 
cleon source operator allows us to explore trial sources with different spatial sizes for the diquark 
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Figure 7: Array of plots of overlap as a function of 
quark displacement, for four-component spinors. 



Figure 8: Array of plots of overlap as a function of 
quark displacement, for two-component spinors. 



and the remaining quark. Contour plots of the resulting two-dimensional parameter space are 
shown in Figs. 5 and 6. The primary feature of these results is the location of the peak overlap 
(very nearly) along the previously explored diagonal cross-section. Since any significant asymme- 
try in the peak would have been suggestive of a strong diquark substructure, this is one piece of 
evidence against the presence of large diquark correlations. 

Relative Displacement of Quark, Diquark 

To explore the physical "dog bone picture" of a separated quark and diquark joined by a flux 
tube, we also tried displacing the single quark from the diquark in our trial sources. Results for 
a variety of displacement lengths and smearing combinations are shown in Figs. 7 and 8. The 
maximum overlap is observed for zero displacement, again suggesting the lack of strong diquark 
correlations. 

4. Conclusions 

In summary, we observe dramatic increases in the overlap between a trial state and the physi- 
cal nucleon as we vary accessible features of the trial state. Smearing the quark interpolating fields 
from a point to an optimal RMS radius increases the overlap for a four-component trial function 
from a fraction of a percent to 35%. Removing the unphysical S-wave lower components increases 
the overlap from 35% to 50%. Optimal smearing of the gauge field further improves the overlap 
from 50% to more than 80% Attempts to further increase the overlap by including diquark corre- 
lations associated with the relative size and position of the quarks and diquarks yielded no further 
improvement, indicating that such correlations are not significantly favored in the nucleon ground 
state. 
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From the perspective of lattice calculation technology, these results are extremely useful in 
generating sources that involve minimal contaminants from excited states. From a physics per- 
spective, they give useful insight into both what the quark and gluon degrees of freedom are doing 
and what they are not doing. It would be valuable to extend these calculations to lighter quark 
masses, but for our research exploiting light domain wall quarks, we need to find an alternative to 
the transfer matrix corrections described in this work. 
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